A Family of Grover's Quantum Searching Algorithms 
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We introduce the concepts of Grover operators and Grover kernels to systematically analyse Grover's 
searching algorithms. Then, we investigate a one-parameter family of quantum searching algorithms 
of Grover's type and we show that the standard Grover's algorithm is a distinguished member of this 
family. We show that all the algorithms of this class solve the searching problem with an efficiency 
of order O(viV), with a coefBcient which is class- dependent. The analysis of this dependence is a 
test of the stability and robustness of the algorithms. We show the stability of this constructions 
under perturbations of the initial conditions and extend them upon a very general class of Grover 
operators. 

PACS number: 03.67.Lx, 03.67.-a 



I. INTRODUCTION 

The problem of searching an element in a list of N 
unsorted elements when this number becomes very large 
is known to be one of the basic problems in Computa- 
tional Science. Classically, one may devise many strate- 
gies to perform that search, but if the elements in the 
list are distributed with equal probability, then we shall 
need to make 0{N) trials in order to have a high confi- 
dence of finding the desired element, also called marked 
element. The formulation of quantum computation as a 
well-stablished theoretical discipline for storing and pro- 
cessing information ||l[ has opened the possibility of de- 
signing new searching algorithms with no classical ana- 
logue. The familiar Grover quantum searching algorithm 
takes advantage of the quantum mechanical properties to 
perform the searching problem with an efficiency of order 

o[^/N) i, i. 

In classical computation there exist design techniques 
which provide general directions for algorithmic problem 
solving. In quantum computation however, the list of 
quantum algorithms is very short. It seems we are lack- 
ing the basic principles underlying the quantum algo- 
rithmic design. Under these circumstances it is a good 
choice to put to the test the currently known quantum 
algorithms. In this work our purpose is to follow this 
goal with Grover's quantum searching algorithm by try- 
ing to understand the relevant pieces of this algorithm 
and wondering to what extent they allow for generaliza- 
tion |, @, |], §, @, §. 

Let us state the searching problem in terms of a list 
£[0, 1, . . . ,N—1] with a number TV of unsorted elements. 
We shall denote by xq the marked element in C that we 
are searching for. The quantum mechanical solution of 
this searching problem goes through the preparation of 
a quantum register in a quantum computer to store the 
N items of our list. This is how quantum parallelism is 
realized. Thus, let us assume that our quantum registers 
are made of n qubits so that the total elements we have 
are N = 2". Let us denote by \x), x — 0,1, . . . ,N ~ 



1 the ket states of the computational basis which are 
orthonormalized. Any state |^) of the quantum register 
is a linear superposition of the computational states. In 
the beginning of the algorithm the quantum register is 
initialized to a given quantum state |^) = |a;in). 

The second component of the algorithm is to design a 
quantum operation which will be repeatedly applied to 
jxin) in order to find the marked element. This strat- 
egy is similar to the classical counterpart algorithm. The 
difference is the fact that the quantum operation is re- 
alized in terms of an unitary operator which implements 
the reversible quantum computation. It is this quantum 
operation what has been so neatly designed by Grover 
0. With Grover's choice we may say that the quan- 
tum evolution is such that the constructive interference 
of quantum amplitudes is directed towards the marked 
state one looks for. 



II. GROVER OPERATORS 

In order to set up our analysis we shall need to intro- 
duce some definitions. 

Definition 1 A Grover operator G is any unitary operator 
with at most two different eigenvalues; i.e., G a linear 
superposition of two orthogonal projectors P and Q: 



G = aP + f3Q, P^=P, Q^ = Q, P 



= 1 



(1) 



where a, (3 € C are complex numbers of unit norm. 
Definition 2 A Grover kernel K is the product of two 
Grover operators: 



K = G2G1 



(2) 



Some elementary properties follow immediately from 
these definitions. 

Property 1 Any Grover kernel i^ is a unitary operator, 
and therefore, it can be used to implement the unitary 
evolution in a quantum computer. 



Property 2 Let the Grover operators Gi,G2 be chosen 
such that 

Gi = aE^g + j3Q^„, P^g = |xo)(a;o|, P^o + Qxo = 1 

(3) 



G2^jP + SQ, P + Q = l 
with P given by the rank 1 matrix 

1 ... 1 



1 



(4) 



(5) 



1 



1 



This is clearly a projector P 
spanned by the state \ko) — 



- |fco)(fco| on the subspace 
^(1,... ,1)*, where the 
superscript denotes the transpose. Then, if we take the 
following set of parameters, 



'1, /3 = 1, 7^ 



-1, (5 = 1 



(6) 



the Grover kernel \^ reproduces the original Grover's 
choice. This property follows inmediately by construc- 
tion. In fact, we have in this case Gi = 1 — IPxa =■ Gx^, 
whilst the operator G2 = 1 — 2P coincides with the dif- 
fusion operator introduced by Grover to implement the 
inversion about the average H. 

One can also show the following property which pro- 
vides a geometrical meaning for the Grover kernels. 
Property 3 Let /Cg denote the set of all the Grover kernels 
for fixed {|a;o), |^o)}- Then /Cq can be viewed as a 3D- 
subset of the group U(2) which is of the form S^ x /Cq, 
where /C^ is a 2D-submanifold of SU(2) (Fig.l). 





FIG. 1. Left: half of the surface /Cq for N = 10. Right: 
the whole surface ¥^q for A^ = 10. In red, the rotation axes 
corresponding to id, 2 £ SU{2). In pale blue, the normal 
to their plane. The yellow dot signals the original Grover 
kernel. The blue curves are the border a = n, Svr, where a 
is the rotation angle. In green, the locus of kernels K with 
factors 16*1.2 with the same rotation angles mod 2-k. 



The content of this property is illustrated in Fig.l. It 
follows from the fact that the two parameters of a Grover 
kernel in Definition 2 with the fixing (|6|) can be used to 
parametrize a subset of the unitary group U(2) of com- 
plex 2x2 matrices. This 3D-subset has a factorized 
form S^ X /Cq where S^ is the unit circle (the group U(l) 



and we call /Cq a certain 2D-submanifold of the group 
of special unitary matrices SU(2) whose construction we 
explain in the following and we plot in Fig.l. This figure 
is constructed by parametrizing the two elements ±iGi,2 
of SU(2) as follows, ±iGi,2 = e'"!.^"!-^'", where a are 
the Pauli matrices and ni.2 are unit vectors which are 
kept fixed. Likewise, we parametrize the corresponding 
Grover kernel (||) as X = e'""'^. Then, upon varying the 
parameters ai,2 we obtain the surface depicted in Fig.l 
Let us point out the following interpretation of the 
Grover operators. Let us think of the computational 
basis {|a;)} as a coordinate basis in Quantum Mechan- 
ics and introduce the quantum discrete Fourier trans- 
form in the standard fashion, la;) :— U-dft\x) = 



^N- 



-m ^y=Q e^^'^s'/^ly). The transformed basis {|x)} can 
then be seen as the dual momentum basis. Then, it is 
easy to see that in such a basis the projector operator P 
takes the following form: 



^DFT^^DFT 



|0>(0|-:Po 



(7) 



This means that the Grover operator G2 takes the same 
matricial form in the momentum basis as the Grover op- 
erator Gi in the coordinate basis. They are somehow 
dual of each other. The original Grover kernel takes then 
the form 



K ~ L/DFTGa;=O^DFT^a: 



(8) 



which shows that a Grover kernel has a part local in coor- 
dinate space and anoter part which is local in momentum 
space. 

This "momentum" interpretation of the search algorithm 
stems from a quantum mechanical analogy between the 
computational basis and its Fourier transformed states 
which enter the definition of the Grover operators. We 
would like to point out that similar analogies have been 
used in connection with alternative formulations of the 
quantum searching algorithm, namely, the analog ana- 
logue of a digital quantum computation with Grover's 
algorithm pC|] which is based in a Hamiltonian formula- 
tion. 



III. THE SEARCHING ALGORITHMS 

A. The Basic Formalism 

Next, the third part of the algorithms corresponds to 
applying the Grover kernel K to the initial state Iccin) a 
number of times m seeking a final state \xf) given by 



|xf(™))=K™|a;in) 



(9) 



such that the probability V{xo) of finding the marked 
state is above a given threshold value. We shall take this 
value to be 1/2, meaning that we choose a probability 



of success of 50% or larger. Thus, we are seeking under 
which circumstances the foUowing condition 



P(xo) = |(xo|i^"|a;in)P>l/2 



(10) 



holds true. 

The analysis of this probability gets simplified if we re- 
alize that the evolution associated to the searching prob- 
lem can be mapped onto a reduced 2D-space spanned by 
the vectors 



{\xq),\x±) 



y/N 



b ^ "^^^ 



X^Xo 



Then we can easily compute the projections of the Grover 
operators Gi , G2 in the reduced basis with the result 



Gi 



a 

/? 



fl2) 



This in turn can be casted into the following closed form: 



(xo|if™|xi, 



1 



^TT + (e" 



1){xo\k2){k.2\Xu 



(18) 



with Aw = W2 — wi. 

In terms of the matrix invariants 



T)eiK = p5, Tri^ = -(/3 + (5) + (l + /3)(l + (5) — 



n 1) the eigenvalues ^1,2 — e'"^'^ are given by 



C1.2 = (Trif )/2 T V-Deti^ + (Tri^/2)2 
The corresponding unnormalized eigenvectors are 



/ A^^-4(Doti<')Af2+A2 
\ki,2) OC 2(1+SWN-1 



(19) 



(20) 



(21) 



Go 



5 

7 



+ (7 - ^) 



1 


\/W-l 


N 


N 


VW-i 


-1 


N 


N 



(13) 



From now on, we shall fix two of the phase parameters 
using the freedom we have to define each Grover factor 
in (g) up to an overall phase. Then we decide to fix them 
as follows: 



a 



7 



-1 



With this choice, the Grover kernel I 
form in this basis: 



(14) 
takes the following 



K 



_ 1 + 5{\ - N) -p{l + S)VI^r^l 
Nyil + 6)yW^^ I3{l + S-N) 



(15) 



We shall fix the initial conditions using the same ini- 
tial state |a;in) as in the original Grover's algorithm 0, 
i.e., we choose the uniform state corresponding to zero 
momentum and find its components in the reduced basis 
to be 



1 



N 



2^0 



N -1. 



N 



(16) 



In order to compute the probability amplitude in (10), 
we introduce the spectral decomposition of the Grover 
kernel K in terms of its eigenvectors {|ki),|k2)}, with 
eigenvalues e'^^ , e"^^ . Thus we have 

A{xo) := (xo|if™|xi„> = 

(17) 



with 



A 



(f3~S)N+il~P)il + S) 



(22) 



Although we could work out all the expressions for a 
generic value N of elements in the list, we shall restrict 
our analysis to the case of a large number of elements, 
N —> 00, and we shall leave for a numerical simulation the 
effect of arbitrary N. Thus, in this asymptotic limit we 
need to know the behaviour for A^ ^ 1 of the eigenvector 
1^2) which turns out to be 



1^2) OC 



13-S 
1+S 



^ + 0(:A,) 



N' 



1 



(23) 



Thus, for generic values of /3, (5 we observe that the first 
component of the eigenvector dominates over the second 
one meaning that asymtoptically \k2) ~ |a;o) and then 
{xQ\K2){K2\xin) = ^(~^)- This implies that the prob- 
ability of success in (Q) will never reach the threshold 
value ([lO|). Then we are forced to tune the values of 
the two parameters in order to have a well-defined and 
nontrivial algorithm and we demand 



(3 = S^-1 



(24) 



Now the asymptotic behaviour of the eigenvector 
changes and is given by a balanced superposition of 
marked and unmarked states, as follows 



^2 



V2\ 1 



(25) 



This is normalized and we see that none of the compo- 
nents dominates. When we insert this expression into 
(H) we find 



\{a:o\Knxin)\-% 



-II 



mALu 



This expression means that we have succeded in finding 
a class of algorithms which are apropriate for solving the 
quantum searching problem. Now we need to find out 
how efficient they are. To do this let us denote by M 
the values of the time step m at which the probability 
becomes maximum; then 



M = Lk/At^lJ 



(27) 




SOO ^ JOOO 



FIG. 2. Probability of success P as a function of the time 
step for N = 1000 and 13 = S ^ e'''/^ 



As it happens, we are interested in the asymptotic be- 
haviour of this optimal periods of time M. From the 
equation ( |20|) we find the following behaviour as iV — )■ oo: 



N 



ReVS 



(28) 



Thus, if we parametrize 6 = e"^, then we finally obtain 
the expression. 



M ' 



IT 



4 cos ■ 



(29) 



Therefore, we conclude that the Grover algorithm of 
the class parametrized by </> is a well-defined quantum 
searching algorithm with an cfhciency of order 0{\/N) 
and with a subdominant behaviour which depends on 
each element of the family. Within this class, the original 
Grover algorithm is a distinguished element for which 
the coefficient in (g9|) achieves its optimal value at = 
0. Moreover, the worst value occurs for — > tt; in this 
limit M is not well-defined, and it corresponds to trivial 
case where the Grover kernel is just the identity operator, 
K = 1. 



The expression (g9|) for M can also be given another 
meaning regarding the stability of the Grover's case 
(26) (/> = 0. It is plain that under a small perturbation 5(f) 
around this value, its optimal nature is not spoiled in 
first order for we find a behaviour which is quadratic in 
the perturbation, namely, M - f (1 + 0.125((5</))2)^/lV. 
This stability considered here is with respect to pertur- 
bations in eigenvalues (or eigenvectors) in the reduced 
2-dimensional subspace specified by the quantum search- 
ing problem (|ll|). We also require these type of pertur- 
bations to hold in all iterations. 

Howewer, if we happen to choose a Grover kernel with a (j) 
far from we may end up with a searching algorithm for 
which the leading behaviour order 0{vN) is masquer- 
aded by the big value of the coefficient and the time to 
achive a succedding probability becomes very large. For 
instance, we may have a Grover kernel with a behaviour 
M ~ 10^-\/iV and for a value of A^ = 10^ it would turn out 
as efficient as a classical algorithm of order 0{N) = 10^. 
Thus, the limit ^ tt behaves as a sort of classical limit 
where the quantum properties disappear. 
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FIG. 3. Probability of success V as & function of the time 
step for iV = 1000 and /? = i, (5 = ie'^^" (up), /? = i, 5 = ie'^ 
(down). 

In Fig. 2 we have plotted the probability of success as 
a function of the time step m for a list of A = 1000 
elements and a choice of parameters (3 = 5 = ^'"1'^ sat- 
isfying condition (24). We observe how the algorithm is 
fully efficient in achieving the maximum probability pos- 
sible. Despite = 7r/2 is not close to Grover's optimal 
value of 0, we find an excellent behaviour. The main 



difference with the optimal case is that here the number 
of maxima is 14 while for Grover's it is 20, as implied 
by (29). This looks like a pattern of fully constructive 



interference. 

In Fig. 3 we have plotted the same function but with 
a choice of parameters /? = i, (5 = ic'^/'* (/3 = i, (5 = ie'"^) 
violating condition (24). We observe how the algorithm 
becomes inefficient and the maximum it takes is lesser 
than 0.0021923 (0.001864) for any time step. This looks 
like a pattern of partially constructive interference. 

We find this behaviour as reminiscent of a quantum 
phase transition where the transition is driven by quan- 
tum fluctuations instead of standard thermal fluctua- 
tions. In this type of transition each quantum phase is 
characterized by a ground state which is different in each 
phase. It is the variation of a coupling constant in the 
Hamiltonian of the quantum many-body problem which 
controls the occurrence of one quantum phase or another 
in the same manner as the temperature does the job in 
thermal transitions. In our case we may consider the two 
different asymptotic behaviours of the eigenvector \k,2) as 
playing the role of two ground states. Following this anal- 
ogy, we may see our family of algorithms parametrized 
by a torus T = S^ x S^ where the parameters /? and 5 
take their values and the difference g :— (3 ~ 5 \s a. sort 
of coupling constant which governs in which of the two 
phases we are. When g 7^ we fall into a sort of disor- 
dered phase where the efficiency of this class of Grover's 
algorithms is spoiled. However, when g = we are lo- 
cated precisely at one equal superposition of the pricipal 
cycles of the torus which defines a one-parameter family 
of efficient algorithms. 



B. The Influence of Initial Conditions 

Next we shall address the issue of to what extent 
this one-parameter family of algorithms depends on the 
choice of initial conditions for the initial state jxin). We 
would like to check that the stable behaviour we have 
found is not disturbed under perturbations of initial con- 
ditions. 

Let us consider a more general initial state |a;in) which 
is not the precise one used in the original Grover's algo- 
rithm PI but instead it is chosen as 



where now |a;in) is the new initial state (30). We have to 
distinguish two cases: i) The coefficient a of the marked 
state is order 1 and ii) it is order bigger than 1, say of 
order 0{^fN). In the latter case ii), it means that the 
initial state is so peaked around the marked state that 
we do not even need to resort to a searching algorithm, 
but instead measure directly on the initial state to find 
sucessfully the marked state. Thus, we shall restrict to 
case i) in the following. Now the key point is to realize 
that all the previous asymptotic analysis is dominated by 
the behaviour of the eigenvector \k2) given by expression 
( P3f ) which is something intrinsic to the Grover kernel 
and independent of the initial conditions. Thus, if con- 
dition ( p4| ) is not satisfied, then as we are in case i) the 
first term in the RHS of (pTl) is not relevant and we are 
led again to the conclusion that the algorithm is not ef- 
ficient. On the contrary, if condition (p3) is satisfied the 
same mechanism based on (G3) operates again and the 
algorithm has a probability of success measured by 



|(xo|X"|xi„)| 



. toAw 
l^|sin(^^— ) 



(32) 



with IS.L0 also given by (28). Then we may conclude that 
the class of algorithms is stable under perturbations of 
the initial conditions. 



C. Extended Formalism 

Finally, we would like to check how general is this con- 
struction in terms of projection operators of the type 
used in (5) for P. To this end let us recall that P can 
be interpreted as the projector |0)(0|. Thus a natural 
generalization is to consider a projector on a different 
momentum state, say \yo), with ya ^ 0. The matrix ele- 
ments of this projection operator in the coordinate basis 
are 



{P)x,.' = le2-(-'--)-!/o/A^^ x,x' =Q,l,...,N^l 

(33) 

We can go even further and consider a general form for 
the states jcco), \x±), |fco) as follows, 



|a;in) = —i=\xo) + o\l — z^\x \ 



N 



(30) 



where a and h are chosen to satisfy a normalization condi- 
tion. Then, it is possible to go over the previous analysis 
and find that the probability amplitude is now given by 



(a;o|i^"Ni, 



^^ + (e' 



lAii 



l)(a;o|K2)(K2|a;in) 



(31) 



|xo) = (l,0,...,0)* 

\xi_) = (0, a2, . . . , unY 

\ko) = (ai,a2, • ■ • ,aAr)* 



(34) 



where there is no loss of generality by chosing |xo) in 
this way; ai, . . . ,aN is a given and normalized set of 
arbitrary complex amplitudes, with ai > 0. We will 
assume that ||ajp > a^. 
The projector P is chosen to be 



P^\ko){ko\ 



(35) 



and it admits (BSh as a particular case. 

Now in the reduced 2D-basis spanned by {\xo), \x±)} 
the Grover kernel has the following expression: 



K = 



-S 



Aaiy^T~. 



-/3Aaiyr 
/3(Aa? - 



1) 



(36) 



with A := 1 + ^. Thus all the dynamics depends on the 
relative strength of the real amplitude ai with respect to 
the rest of the amplitudes. If we set at = 1/VN Vi, then 
we recover the same expression as in dlq). Moreover, the 
initial condition is taken as 



Ixin) = ai|a;o) + A/l - al\x±) 



(37) 



In order to perform our analysis, we shall assume that 
the unknown amplitude ai behave s generically as ai 



1/VN, and consequently i/l — af ~ ^yl — 1/N. Under 
these circumstances, we find the following asymptotic be- 
haviour for the eigenvector |k2) of the Grover kernel: if 
(3^5 and 6^ -1, 



(38) 







/ /3-<5 1 
K2) OC ( 1+^ "1 


if/3 = 


= -5, 


1 X 1 ^i<5'/' 



(39) 



This latter case is again the only favorable to obtain an 
efficient algorithm and the behaviour of the time M for 
achieving maximum probability of success takes the fol- 
lowing form 



M ^ 



Tran 



4 cos- 



(40) 



We conclude then that our construction of quantum 
searching algorithms of Grover's type are general enough 
under different choices of Grover operators Gi , G2 and 
that the analyisis performed with the simplest choice of 
these operators captures the essential properties of the 
class of algorithms we have presented. 



facilitates the generalization of Grover's algorithms in 
several direcctions. We have characterized the basic fea- 
tures of these algorithms in terms of these operarators 
whose main properties we have established in Sect. II. 
Using these operators we have investigated a family of 
Grover kernels whose qualities as efhcient algorithms de- 
pend on the range of parameters entering the construc- 
tion of their associated Grover operators. When the al- 
gorithms are efficient, they also perform the searching 
task with order 0{^/N), and the original Grover's choice 
gives the optimum value in the one-paramater family of 
algorithms. Moreover, we have extended this study to 
incorporate initial conditions different than the standard 
uniform initial states and we have checked that letting 
aside exceptional cases, the basic algorithms of Sect. Ill 
maintain their efficiency. Finally, we have addressed also 
the issue of considering quite general Grover operators 
and found that the basic efficiency properties of the sim- 
plest choice's for Grover's algorithm remain unchanged. 
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